Collective modes in relativistic npe matter at finite temperature 
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Isospin and density waves in neutral neutron-proton-electron (npe) matter are studied within a 
relativistic mean-field hadron model at finite temperature with the inclusion of the electromagnetic 
field. The dispersion relation is calculated and the collective modes are obtained. The unstable 
modes are discussed and the spinodals, which separate the stable from the unstable regions, are 
shown for different values of the momentum transfer at various temperatures. The critical tem- 
peratures are compared with the ones obtained in a system without electrons. The largest critical 
temperature, 12.39 MeV, occurs for a proton fraction j/p = 0.47. For yp = 0.3 we get Tcr=5 MeV 
and for yp > 0.495 Tcr S; 8 MeV. It is shown that at finite temperature the distillation effect in 
asymmetric matter is not so efficient and that electron effects are particularly important for small 
momentum transfers. 
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I. INTRODUCTION 

In order to describe compact-star matter at low den- 
sities, besides protons and neutrons, the electrons (and 
possibly neutrinos if they are trapped) must be intro- 
duced. Electrons neutralize the proton charge and thus 
suppress the diverging Coulomb contribution to the en- 
ergy. This composite matter is expected to undergo 
phase transitions associated to the nuclear liquid-gas 
phase transition but strongly modified by the Coulomb 
interaction and the presence of electrons. Phase tran- 
sitions in asymmetric matter are related with the in- 
stability regions limited by the spinodal sections 0, Q] ■ 
These very same instabilities are responsible for nucle- 
ation processes Q and the possible existence of a non- 
homogeneous phase. At low nuclear matter densities, 
which are the densities of relevance in the present work, 
a competition between the long-range Coulomb repulsion 
and short-range nuclear attraction can lead to the forma- 
tion of matter known as nuclear pasta Q. The compe- 
tition of these two forces, normally acting on different 
scales, is called frustration. This novel state of matter, 
the nuclear pasta, can appear in different structures and 
its properties have important consequences in the crust 
of neutron stars and in the core-collapse of supernova 
0, H, Q, I3- Not only the equation of state of stellar 
matter has to be understood, but also the neutrino mean 
free path in the medium has to be well described. In fact, 
neutrino interactions are crucial in the dynamics of the 
core-collapse supernovae because they carry most of the 
energy away. It has been shown that the neutrino opacity 
is affected by nucleon-nucleon interactions due to coher- 
ent scattering off density fluctuations 0- Both single 
particle and collective contribution have to be taken into 



account. In [lOj it was shown that coherent neutrino scat- 
tering from non-uniform hadron-quark matter or hadron 
matter with and without kaon condensed phase would 
greatly reduce the neutrino mean-free path. A similar 
effect at low densities could allow enough energy transfer 
in order to revive the supernova shock. Recent semiclas- 
sical simulations of the linear response of nuclear non- 
homogenous matter at low densities, the so called nuclear 
pasta, to neutrinos in (51, Q have shown that coherence 
effects reduce the mean free path of neutrinos. In these 
simulations electrons are not modeled explicitly, but their 
effect is included through a modified Coulomb interaction 
between the protons through a screening length. Neutri- 
nos will loose energy by exciting collective nuclear modes 
or plasmon modes. However in it has been shown 
that the behavior of the electrons depends on the wave- 
length of the perturbation. In low energy nuclear 
collective excitations of Wigner-Seitz cells containing nu- 
clear clusters immersed in a gas of neutrons have been 
obtained. In that paper, the electron motion was not 
included. Including the electron contribution will likely 
affect the results for large clusters. 

In ^3 we have investigated the influence of the elec- 
tromagnetic interaction and the presence of electrons on 
the unstable modes of nuclear matter and compared the 
dynamical instability region with the thermodynamical 
one. This investigation was performed in the framework 
of a relativistic mean field hadronic model within the 
Vlasov formalism 

0, E El- In particular, in pdj. it 
has been studied the role of isospin and the modifica- 
tion of the distillation phenomenon due to the presence 
of the Coulomb field and electrons. These calculations 
were performed at T = MeV. 

In the present work the effect of dynamical instabilities 
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including the Coulomb field are investigated at finite T. 
Understanding the mixed phase of neutral matter at fi- 
nite temperature is important to determine the behavior 
of neutrinos emitted in a supernova explosion. 



II. THE VLASOV EQUATION FORMALISM 

We consider a system of baryons, with mass M inter- 
acting with and through a isoscalar-scalar field cj) with 
mass TOs, a isoscalar- vector field with mass m„ and 
an isovector- vector field b'^ with mass irip. We also in- 
clude a system of electrons with mass me. Protons and 
electrons interact through the electromagnetic field A^^ . 
The lagrangian density reads: 

C=Y^ C^+Ce+Ca+C^+Cp+CA. (1) 

where the nucleon Lagrangian reads 

A =i^,[7m^^^-M*]V,, (2) 

with 



^ 2 2 

M* = M-gs^, 

and the electron Lagrangian is given by 

= 4 [if. + eA^) - me] Ve- 



(3) 
(4) 

(5) 



The isoscalar part is associated with the scalar sigma (a) 
field, (p, and the vector omega (w) field, Vp, while the 
isospin dependence is coming from the isovector vector 
rho (p) field, (where /i is the 4 dimensional space- 
time indices and i the 3D isospin direction indices). The 
associated Lagrangians are 



= 

r — 

= 

Ca = 



-^n^M'^-' + \mlVpV" 
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where Vt^i, = d^Vy - d^Vp , B^y = d^^by - dybp - gp{bp x 
by) and Fp_y — dpAy — dyAp. The model comprises the 
following parameters: three coupling constants gs, gv and 
gp of the mesons to the nucleons, the nucleon mass M, 
the electron mass me, the masses of the mesons m^, m^, 
nip, the electromagnetic coupling constant e = y/ 47r/137 
and the self-interacting coupling constants k and A. We 
have used the set of constants identified as NL3 taken 
from For this case, the saturation density that we 



refer as po is 0.148 fm ^. In the above lagrangian density 
T is the isospin operator. 
Denoting by 

/(i",P,i)± = diag(/p±, fn±, fe±) 

the distribution functions of particles at position r, 
instant t and momentum p and of antiparticlcs (-) at 
position r, instant t and momentum — p and the corre- 
sponding one-body hamiltonian by 



h± = diag (hp±,hn±,he±) 

where 

For protons and neutrons, i — p,n, we have 
Voi = gvVo + y Tibo + eAo^-^, 



V.=g,„V+|r,b + eAi±^ 



(6) 



M; = M* = M* = M - g,q}, 

Ti = 1 (protons) or -1 (neutrons), and for electrons, i = e, 
we have 

Voe = -eAo, Ve = -e A, M* = m^. 

The time evolution of the distribution functions is de- 
scribed by the Vlasov equation 



i± 



dt 



I = p, n, e. 



(7) 



where {, } denote the Poisson brackets. From Hamilton's 
equations we derive the equations describing the time 
evolution of the fields (f>, V", A^^ and the third compo- 
nent of the /9-field 63 = (69, b), which are given in the 
Appendix A. 

The state which minimizes the energy of asymmetric 
npe matter is characterized by the distribution functions 



1 



I — p,n 



with e = y/jP+M* 



t^i - 9vVo - y Tibo - e Aq — - — 



and 



1 



ie± 



1 -|_ g(eoeTPe)/T 



(8) 



(9) 



(10) 
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with eoe = \/p^ + TOg, and by the constant mesonic fields 
which obey the following equations 



2 T^(0) .(0) 



^2,(0) _ 5p .(0) 



(0) 



l(0) 



A 



(0) 



0. 



Collective modes in the present approach correspond 
to small oscillations around the equilibrium state. These 
small deviations are described by the linearized equations 
of motion and, therefore, collective modes are given as 
solutions of those equations. To construct them, let us 
define: 

fi± = foi± + <5/i± , 



bo = 60"'' + Sbo , bi = Shi 



i = p,n, where 



SVoi = gvSVo + Sbo + e '^^0, 



SV^ = gvSV + n^6h 



e — - — oA, 



with hoi± = ±e + Vq, , i = p,n and hoe± — ±eoe ■ 
The linearized equations of the fields are shown in the 
Appendix B. 

Of particular interest on account of their physical rele- 
vance are the longitudinal modes, with momentum k and 
frequency w, described by the ansatz 




V 



5BI 



where j = e, p, n, B = V, 6, A represents the vector fields 
and 9 is the angle between p and k. For these modes, we 
get 5V^ = Syy ^0,6bZ=Sby^0 and 6AZ = = . 
Calhng SV^ = SV^, Sb^ = Sb^ and SA^ = SA^, we will 
have 

and the equations of motion become 

1 rTiQ f) 

, (13) 



kpcosO 



Ao = (5Ao , A^ = (5ylj 



(11) 



As in milTl we express the fluctuations of the distri- 
bution functions in terms of the generating functions 



kp cos 9\ : Af* (. , pcos6'-, „■ ^, ,n,■ 

(14) 



S'±(r, p, t) = diag (Spi, S'„±, S'e±) , 



such that 



'^/± = {5±,/o±} = {5±,p^} 



4fo=i 



In terms of the generating functions, the linearized 
Vlasov equations for 5]i± 



d6fi± 
dt 



+ {Sfi±, ho,±} + {foi±,Shi±} = 



are equivalent to the following time evolution equations 



e± 



dt 



{Se±,hoe±} — She± — — £ 



MoT 



P-6A 

eoe 



— — V{bi±,hoi±} = dhi± = =F 

at e 



1 

\l2) 



_ fc2 _ ^2^^^^) 



AigskM* ^ 5ps^ (15) 



{oj^ ~k^~ ml) 5V° = 4ig,k > (Sp, (16) 



E 



■^-e-ml)5bl = AY,TM (17) 



fc2)M° = 4iefc ^(-l)"'<5p„ (18) 



where, np = 0, ng = 1, and, from the continuity equation 
for the density currents, we get for the components of the 
vector fields 



(19) 
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(20) 



IV. NUMERICAL RESULTS AND 
DISCUSSIONS 



In the above equations we have defined 
with 



(21) 



(22) 







Sps 



(27r)3 

d^p pcosO 
(27r)3 e 



(/oi+ + /w-)^ 



5, 



dp^ 
dfoi- 



dp"^ 



III. 



^ - ^/o,±(/o.± - 1). 



SOLUTIONS FOR THE EIGENMODES AND 
THE DISPERSION RELATION 



The solutions of eqs. — H21I) form a complete set 
of eigenmodes which may be used to construct a general 
solution for an arbitrary longitudinal perturbation. 

Substituting the set of equations (|15I21|) into (|13|l and 
H14|l we get a set of equations for the unknowns Sl^^ 
which lead to the following matricial system 



/ an ai2 ai3 au ai5 \ 

021 022 023 024 

031 032 033 a34 a35 

041 042 043 a44 

V a53 / 



Pujp 

PuJ7l 

\ Pu>e / 



0. 



The coefficients defined in the Appendix C. The 

amplitudes j and p^i are functions of the quantities 
S'ii given in Appendix C, Eqs. (IHIJ and The 
dispersion relation is obtained from the determinant of 
the matrix of the coefficients, namely 

On 012 oi3 ai4 ai5 
021 022 023 024 
O31 032 033 034 035 

O41 042 043 O44 
053 O55 

In the next section we discuss the behavior of the ratios 
of amplitudes Pt^n/ puip and Pue/ Puip at or inside the spin- 
odal region as a function of the density and the isospin 
asymmetry. 
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FIG. 1: Spinodal sections for several temperatures and mo- 
mentum transfer fc = 75 MeV and fc = 20 MeV. 

In this section we present the results for the modes that 
define the instabilities of the system, corresponding to 
imaginary frequencies. We are interested in understand- 
ing the effect of the temperature and of the presence of 
electrons on these instabilities. 

In Fig. ^we plot the spinodal region for different tem- 
peratures and k = 20, 75 MeV. The second value of k 
approximately defines the envelope of the spinodal re- 
gion. We have also considered a smaller value of k in 
order to get a larger effect of the electrons. As can be 
seen in Fig. O the influence of the electrons decreases, as 
k increases up to /c < 70 — 80 MeV/c (thick lines), similar 
to what happened at zero temperature. This is expected 
since the Coulomb contribution varies with the inverse of 
the momentum square, becoming weaker at large k. For 
higher values of /c, this effect becomes negli gibl e and we 
recover the behavior already observed in 0,[ill, i.e., the 
instability region decreases with the increase in the mo- 
mentum transfer (thin lines). This is in fact coming from 
the finite range of the nuclear interaction which reduces 
the binding of the matter with a fc^ term and leads to 
a reduction of the spinodal region at large k. In a ther- 
modynamical calculation of the spinodal for npe matter 
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FIG. 2: The ratio |^ at the critical points. Thin (thick) hnes 
represent the critical points at the section pp > p„ (pp < Pn). 



within the NL3 parametrization done in 16] it was shown 
that above 2.4 MeV the instability region completely dis- 
appeared. In fact it is important to include electrons in a 
dynamical way and allow independent electron and pro- 
ton density fluctuations in order to get the real behavior 
of the system. For fc = 75 MeV it is seen that an instabil- 
ity region exists for quite high temperatures, T^r = 12.39 
MeV. The largest critical temperature is a function of 
k and, as we discuss next, occurs for symmetric matter 
only for matter without electrons or for very large values 
of fc, when the Coulomb field is negligible. The spin- 
odals with A; = 20 MeV present several differences with 
respect to the fc = 75 MeV case: the critical tempera- 
ture is lower, the range of densities corresponding to the 
instability region is smaller and, while for fc = 75 MeV 
the spinodals are practically symmetric with respect to 
equal proton/neutron concentrations, for fc = 20 MeV 
the spinodal sections have a larger contribution in the 
Pp < pn range. 

In Fig. n we also include the line of critical points. 
These points of the spinodal are such that the associ- 
ated eigenvector of the equations of motion is tangent to 
the spinodal curve. A very similar method was obtained 
to determine the critical points of the thermodynamical 
spinodal in Again, it is clear the difference between 
the two values of fc discussed. In Fig. |21 we plot the 
ratio 1^ at which the critical point occurs. The lower 
branches (thin curves) for both cases correspond to the 
Pp > Pn section while the upper branches (thick curves) 
to the pp < Pn section. For fc = 75 MeV the ratio of pro- 
ton/neutron density fluctuations is practically parallel to 
the line of equal concentration of protons and neutrons 
(pn/ Pp = l)i not much different from the behavior of np 
matter jl3| ■ However for fc = 20 MeV the eigenvectors at 
the instability border are tilted in the neutron direction, 
preventing a strong distillation effect. 

As already referred, in Fig. 3 we show the spinodal 
regions for different values of fc at T = and 10 MeV. 
We also include the equation of state (EoS) of stellar 
matter in /3-equilibrium without neutrinos {Y^,^ = 0) and 
with trapped neutrinos (Yl 7^ 0). In this last case we 
have considered a constant leptonic fraction (electrons 



plus neutrinos) equal to Yl = 0.4, according to [l^ . 

For T= 10 MeV, an unstable region only exists for 
fc > 30 MeV and fc < 150 MeV. The spinodal region is 
strongly asymmetric at small fc because of the electron 
stabilizing effect, it becomes almost symmetric at large 
fc with the reduction of the coupling with electrons. For 
T= MeV it is also seen the strong asymmetry of the 
spinodal at k=5 MeV as compared to k=250 MeV when 
the effect of electrons is negligible. 

The EoS of stellar matter in /? equilibrium crosses all 
the spinodal sections above fc=6 MeV and below fc=150 
MeV at T=0 MeV indicating the existence of a nucle- 
ation phase at the crust of the star. The size of the 
inhomogeneities is of the order 8-200 fm. However, from 
Fig. |Slwe see that the most unstable mode has fc — 145 
MeV which corresponds to a modulation with A ~ 9 fm. 
At r=10 MeV the /3 equilibrium EoS does not cross the 
unstable region. At this temperature matter is homoge- 
neous for all densities. However, for matter with trapped 
neutrinos the inhomogeneities are again present. In this 
case, the EoS crosses the unstable region even at T=10 
MeV for 30 MeV< fc < 137 MeV and at T^O MeV the 
spinodals are crossed for all values of fc for which the 
instabilities occur. 



T= 10 MeV 
P equilibrium 
Y„ =0 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 



FIG. 3: Spinodal sections at T=0 and 10 MeV. The EoS for 
/3- equilibrium matter is also shown. 
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In Fig. the ratio of the neutron-density fluctu- 
ation over the proton-density fluctuation is plotted as 
a function of momentum for T = and 5 MeV and 
Up — 0.5, 0.3, 0.2 and p — O.lSpo- Comparing the finite 
temperature results with the zero temperature ones we 
conclude: a) for large k the fluctuations are closer to the 
ratio Pn/Pp and therefore, the distillation effect is not so 
efficient: b) the anti-distillation effects, 4^ > — , settles 
in at larger values of k. This becomes more clear from 
Fig. ^p) and it is more pronounced for asymmetric mat- 
ter. In particular for yp — 0.3 the ratio Sp„/Spp comes 
close to the value Pn/Pp = 2.33 (denoted by the dotted 
line) for k — 25 MeV and is completely above that value 
for fc = 15 MeV. It is also interesting to notice that for 
fc = 15 MeV, contrary to fc = 25, 75 MeV, asymmetric 
matter {yp = 0.3) is more unstable at higher temper- 
atures than symmetric matter. This point is discussed 
below together with Fig. |H| 

The most unstable mode (the one with the largest 
modulus of the imaginary frequency), indicating the exis- 
tence of an instability, defines the behavior of the system, 
namely, the most probable size of the inhomogeneities 
which are formed due to the perturbation. In Fig. [51 
we plot w^Q^ and the momentum fc at which this mode 
occurs as a functions of temperature and isospin asymme- 
try. The modulus of w^^j. decreases with temperature, 
isospin asymmetry and the presence of electrons. This 
behavior is expected since any of these factors decreases 
the instability region of the system. 

In Fig. ISJd) we plot k{uJmax) for the values of Fig. [SJi) 
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FIG. 4: The ratio of the neutron-density fluctuation over the 
proton-density fluctuation is plotted for p = 0.15po as a func- 
tion of a) momentum for T — (thick) and 5 MeV (thin); b) 
temperature for fc = 15 (very thick line), 25 (thick line) and 
75 MeV (thin line). In this figure the dotted line represents 
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FIG. 5: Most unstable mode as a function of temperature for: 
a) the square of the corresponding frequency at p = 0.3po 
and yp — 0.2 (black) and yp — 0.5 (grey); b) the momentum 
transfer at which the most unstable mode occurs for several 
densities and proton fractions (same colors as in a) ). 



(p ~ 0.3po and yp — 0.5, 0.2) and for yp — 0.2, p — O.lpo 
and 0.5po- For the same density, the value of k{uj„iax) 
also decreases both with the increase of temperature and 
isospin asymmetry. This means that the effect of the 
electrons, which is more important for smaller fc values, 
is more effective at higher temperature values and larger 
asymmetries. For a given temperature, electrons shift the 
most unstable mode to slightly larger values of fc. If we 
fix the proton fraction, yp = 0.2, and vary the density we 
conclude that close to the critical temperature the fc lies 
always in the interval 60-80 MeV. The large difference 
occurs at zero temperature; in this case the smaller the 
density the larger the value of fc. These results seem to 
indicate that only under special conditions the electrons 
have a strong effect on the behavior of the system, since 
the most unstable mode has fc > 60 MeV and larger 
electron effects are expected for fc below this value as 
discussed next. 

In Fig. iniwe plot the critical temperatures and densi- 
ties as a function of the wave number of the perturbation 
for npe matter and nuclear (np) matter. For values of 
fc > 150 MeV the critical temperatures almost coincide. 
As expected, electrons decrease drastically the critical 
temperatures for small values of fc. Another effect is a 
reduction of the critical density at low fc transfer which 
is more pronounced on the proton density than on the 
neutron density. A high proton density region will have 
a larger contribution from electrons which stabilize mat- 
ter. We conclude that while at np matter the critical 
temperature occurs for symmetric matter, this is gener- 
ally not anymore true if electrons are present. In fact, 
for large values of fc (fc > 150 MeV) matter is almost 
symmetric at T = Tcr, but for small values of fc we can 
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have significant asymmetries for yp ~ 0.25 — 0.4 at the 
critical temperatures. 



V. CONCLUSIONS 

In the present work, we have looked at the dynamical 
instabilities of nuclear neutral warm matter which is of 
interest for the study of neutron stars and supernovae. In 
particular, we have studied the effect of the temperature 
and the inclusion of the Coulomb interaction on the in- 
stability region of nuclear matter. The calculations were 
performed within a relativistic mean-field approach to 
nuclear matter, namely the NL3 parametrization of the 
NLWM but we believe that the main conclusions with 
respect to the dynamical instabilities would be similar 
if obtained within other models with constant couplings 
taking into account the results of . 

We may interpret the present calculation as a way of 
going beyond a total screening effect calculation which 
considers equal distributions of electron and proton 
clouds. On the other hand, no electron screening would 
correspond to an homogeneous distribution of electrons. 
We have seen that for small k transfers, electrons and pro- 
tons stick together which corresponds to a large screening 
effect. On the other hand for large momentum transfers 
protons and electrons behave independently and this cor- 
responds to the no screening limit. 

For np matter it was shown that the largest critical 
temperatures above which nuclear matter is stable oc- 
curs for symmetric matter. This is not anymore true for 
npe matter: critical temperatures occur for asymmetric 
matter with yp < 0.5. The critical proton fraction can 
be as low as 0.22 for fc = 8 MeV, where Tcr ^ 4.5 MeV. 
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FIG. 6: Critical a) temperatures, b) proton and neutron den- 
sities, as a function of momentum transfer including/not in- 
cluding electrons. 



In [g the structure of hot dense matter at subnuclear 
densities has been investigated by quantum molecular 
dynamics (QMD) simulations. In this work the criti- 
cal temperature for the phase separation is >6 MeV for 
the proton fraction yp — 0.5 and >5 MeV for j/p — 0.3. 
These results seem to be compatible with Fig 1^1 In our 
approach we get for yp = 0.3 a critical temperature of 5 
MeV at fc = 10 MeV. The maximum critical temperature 
is 12.38 MeV and occurs for yp = 0.47 at fc = 75 MeV. 
For npe matter with yp > 0.495 the critical temperature 
occurs at T < 8.5 MeV for k > 180 MeV. 

We have seen that stellar matter with trapped neu- 
trinos is more affected by the presence of a non- 
homogeneous phase than neutrino free stellar matter. We 
have also concluded that at finite temperature the dis- 
tillation effect in asymmetric matter is not so efficient, 
and in particular the anti-distillation effect defined as 
4^ > which occurs at small fc, settles in at larger 
values of k. We may expect this effect to have conse- 
quences in warm stellar matter. In a supernovae explo- 
sion 99% of the energy of the collapse is radiated away 
in neutrinos. Neutrinos interact strongly with neutrons 
because of the large weak vector charge of the neutron. 
Therefore, the way neutrons clusterize is important to 
determine the neutrino mean free path. These neutri- 
nos may couple strongly to the neutron-rich matter low- 
energy modes present in this explosive environment |l2j| , 
and in this way revive the stalled supernovae shock. 

The wavelength of the modulations which correspond 
to the most unstable modes increases with temperature. 
While for T = MeV the presence of electrons has al- 
most no effect due to the high values of k (~ 125 — 145 
MeV), for asymmetric matter close to the critical tem- 
peratures they really make some difference. These are 
temperatures and proton fractions of the order of the 
ones found in stellar matter of protoneutron stars and 
therefore will affect the crust formation and the neutrino 
diffusion during the cooling process. 

The results we have obtained are intimately related 
with the studies developed in 0, S 111 j where the nuclear 
pasta is generally modeled as a charge neutral system of 
protons, neutrons and electrons at subnuclear densities. 
In one sees that, for small momentum transfers, the 
dynamical response function shows a peak, identified as 
a plasmon mode, and displays an important amount of 
strength at low energies identified as a coherent density 
wave. In this work a fixed screening length of 10 fm 
was considered. Small momentum transfers are precisely 
the ones for which protons are more influenced by the 
presence of electrons and therefore a carefuU calculation 
taking into account this point should be carried out. 

In ^ the authors have compared different treatments 
of Coulomb interaction: in particular they have shown 
that the full and the no electron screening calculations 
give similar results demonstrating the weakness of the 
electron-screening effects. We came also to a similar con- 
clusion, however we could identify under which condi- 
tions the electron-screening effects are important, namely 
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for small k transfers. 

In it was shown that density dependent hadronic 
models have different behaviors from the NLWM we 
have studied in the present work, coming closer to non- 
relativisitc nuclear models with Skyrme forces. The dy- 
namical instabilities of npe matter within these models is 
currently under investigation. The calculation of collec- 
tive modes of neutron-rich nuclear clusters immersed in a 
gas of neutrons, including electrons dynamically, should 
also be studied. 



APPENDIX A: EQUATIONS FOR THE FIELDS 



joe(r,t) = 2 J ^(/e+(r,p,i)-/e-(r,p,<)) , (All) 



f (fp 



i) + /e-(r,p,t)) 



p -I- eA 



(A12) 

and the components of the isovector four-current density 
are 



J3,0 



i—p,n 

(A13) 



_ V2<^ ml(b + '^<f + = 5sPs(r, t) , (Al) 



(Pp p-Vi 



T, (/.+ (r,p,t)+/,_(r,p,i)) , 
(A14) 



^0 V72 



VVo + m^Vb gvjo{r,t) , (A2) with e, ^ {p - Vi)'^ + M*^ ,i = p,n and = 

^{p + eA)2 + ml . 



VV.-fm^K = ff.Ji(r,t) , (A3) 



APPENDIX B: LINEARIZED EQUATIONS FOR 
THE FIELDS 



^-V^6o + rn>o.|,3,o(r,t), (A4) ^ 

- \7^6(b + [ml + K<j)o + = .9, Sp^ , (Bl) 

^-V\ + m^6. = |j3.(r,i), (A5) 

~ V^SVo + mlSVo = g-a Sjo , (B2) 

^ - V^Ao - e[iop(r, t) - joe(r, t)] , (A6) 

- V^6V, + mlSV, = 6j, , (B3) 

^ - V^A, = e[jp,(r, i) - U{r, t)] , (A7) 2 

y'Sbo + mlSbo = ^ <5j3.o , (B4) 



where the scalar density is 



Jo(r,t) = 2^ J ^(/,;+(r,p,t)-/,_(r,p,t)) , 



Psir,t) = 2j: /^(/.+ (r,p,t) + /._(r,p,t))— . ^^^2,.^^.,. ^ 9j> 



(27r)3-'— ■ ___i_v^56.-fm>. = ^5j3., (B5) 

(A8) 

The components of the baryonic four-current density 
are d 5A^ 2 



- V^JAo = e [Jjop - 5jOe] , (B6) 



(A9) - VHA, = e [Jjj,. - 5j,i\ , (B7) 



= / (|;^(/.(r,P,t) + /._(r,p,t))P-^, / d^P M\,^ , ^ „ 

-P."-" ^ ' ^^^^^ Sps - 2 ^ y ^ —{Sf,+ + 5U^) - gsS<t> dpi , 



i=p,n 



i—p.n 



i=p,n 



i=p,n 



(27r)3 ' e 



eoe 



+ 2e 



(27r) 



3 (/0e+ + /Oe-) ; P^^ 

eOe Coe 



APPENDIX C: DISPERSION RELATION 

The dispersion relation is obtained from the determi- 
nant of the following set of five equations 



\ 






Pion 




pLop 




Pojn 


/ 





0-21 0-21 <I23 <I24 

^31 ^32 ^33 CI34 ^35 
a4i a42 <243 044 



The coeficients a^- are given by 

an = l + c«(/°^+-C_), 



0. 



041 



a2i=c,(/0;-e), 
a22 = l + c« /O'i), 
a23 = -(c.-c,) + 
a24 = -(c.+c^) (4; + ^-), 

(125 = 

a3i = +c. + e) , 

a33 = 1 - (C. + + Ce) (/^^ - /^^_) 
034 = -(C„ - Cp) (7^^ - /^^) 
a35 = Ce - e) 



_ „ / rln I t\ 

O41 = + 

a42 = (7i![ + /il) , 
a43 = -(c.-cp) 
a44 = l-(c.+c^) (/^!^-/^!^), 

(I5I = fl52 = (I54 = 0, 

055 = 1 - Ce - 

where we have defined the following quantities 

2G\ 1 



(27r)2r 0)2 - a;2 
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Cij 



1_ (I)' ^2 ^ (fc2+^2)/fc2^ where 



(27r)2rw2 



(27r) 
pcos^ 



/p/e 



X 



OJ ±x 



± 



2 — \- to In 
e 



u-p/e 



-p/t 



i-oo i-p/e 

AZi^^ = / e"de / dx X {x, p) foi± {foi± - 1) • 

(C3) 

At low densities, corresponding to a negative valnc of 
the compressibiHty, the system presents unstable modes 
characterized by an imaginary frequency In order to 
obtain these modes, one has to replace u) by i/3. In this 
case, the integrals I^j^ (e) become 

l0^{e) = ± [2^ - 2/?tan-i(p/e/3).l (C4) 



The amplitudes and p^^i are given by 

Pa,i = ^°,,+ +^°^,- (C2) 
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